Isotope effect on the electron band structure of doped insulators 
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Applying a continuous-time quantum Monte-Carlo algorithm we calculate the exact coherent band 
dispersion and the density of states of a two dimensional lattice polaron in the region of parameters 
where any approximation might fail. We find an isotope effect on the band structure, which is 
different for different wave- vectors of the Brillouin zone and depends on the radius and strength of 
the electron-phonon interaction. An isotope effect on the electron spectral function is also discussed. 

PACS numbers: PACS: 74.72.-h, 74.20.Mn, 74.20.Rp, 74.25. Dw 



The isotope substitution, where an ion mass M is var- 
ied without any change of electronic configurations, is 
a powerful tool in testing the origin of electron correla- 
tions in solids. In particular, a finite value of the iso- 
tope exponent a = —d\nTc/d\nM proved that the 
superconducting phase transition at T = T^. is driven 
by the electron-phonon (e-ph) interaction in conventional 
low-temperature superconductors. Advances in the fab- 
rication of isotope substituted samples made it possible 
to measure a sizable isotope effect also in many high- 
temperature superconductors. This led to a general con- 
clusion that phonons are relevant for high Tc- Essential 
features of the isotope effect on Tc, in particular its large 
values in low cuprates, an overall trend to decrease as 
Tc increases, and a small or even negative a in some high 
Tc cuprates, could be understood in the framework of the 
bipolaron theory of high-temperature superconductivity 
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The most compelling evidence for (bi)polaronic carries 
in novel superconductors was provided by the discovery 
of a substantial isotope effect on the (super)carrier mass 
0, 3 . The effect was observed by measuring the mag- 
netic field penetration depth Xh of isotope-substituted 
copper oxides. The carrier density is unchanged with 
the isotope substitution of O^^ by O^^, so that the iso- 
tope effect on Xh measures directly the isotope effect 
on the carrier mass to*. A carrier mass isotope expo- 
nent am — dhim* /dlnM was observed, as predicted by 
the bipolaron theory 3- ordinary metals, where the 
Migdal adiabatic approximation is believed to be valid, 
am = is expected. However, when the e-ph interaction 
is sufficiently strong and electrons form polarons (quasi- 
particles dressed by lattice distortions), their effective 
mass TO* depends on M as to* = mexp^jEp/uj). Here to 
is the band mass in the absence of the electron-phonon 
interaction, Ep is the polaron binding energy (polaron 
level shift), 7 is a numerical constant less than 1 that 
depends on the radius of the electron-phonon interac- 
tion, and w is a characteristic phonon frequency (we use 
h = 1). In the expression for to*, only the phonon fre- 



quency depends on the ion mass. Thus there is a large 
isotope effect on the carrier mass in (bi)polaronic con- 
ductors, am — (1/2) ln(TO*/TO) 01, in contrast with the 
zero isotope effect in ordinary metals. 

Recent high resolution angle resolved photoemission 
spectroscopy (ARPES) [6^ provided another compelling 
evidence for a strong e-ph interaction in the cuprates. 
It revealed a fine phonon structure in the electron self- 
energy of the underdoped La2_2;Sr2;Cu04 samples 0, 
and a complicated isotope effect on the electron spectral 
function in Bi2212 P|. Polaronic carriers were also ob- 
served in colossal magneto-resistance manganites includ- 
ing their low-temperature ferromagnetic phase, where an 
isotope effect on the residual resistivity was measured 
and explained 0. 

These and many other experimental and theoretical 
observations point towards unusual e-ph interactions in 
complex oxides, which remain to be quantitatively ad- 
dressed. While the many-particle e-ph problem has been 
solved in the weak-coupling, A = E„/zt <C 1 [sl ITolllTI|. 
and in the strong-coupling A ^ 1 |12j limits, any analyt- 
ical or even semi-analytical approximation (i.e. dynamic 
mean-field approach in finite dimensions) is refutable in 
the relevant intermediate region of the coupling strength, 
A ~ 1, and of the adiabatic ratio, us/t ~ 1. Here t and z 
are the nearest neighbor hopping integral and the coor- 
dination number of the rigid lattice, respectively. 

Advanced variational 0|, direct diagonalization [l^ . 
and quantum Monte Carlo (QMC) [ll IH 113 tech- 
niques addressed the problem in the intermediate region 
of parameters, but mainly in the framework of the Hol- 
stein model ^] with a local (short-range) e-ph interac- 
tion, or in the continuous (effective mass) approximation 
for the bare electron band p^ . However, in cuprates 
there is virtually no screening of c— axis polarised opti- 
cal phonons because an upper limit for the out-of-plane 
plasmon frequency (< 200 cm~^[l9|) is well below the 
characteristic phonon frequency w ~ 1000 cm~^. Hence, 
the unscreened long-range Frohlich interaction is impor- 
tant in cuprates and other ionic charge-transfer insulators 
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El mil mill. A parameter-free estimate ^ yields 
i?p ~ 1 eV with this interaction alone, which is larger 
than a magnetic (i.e. exchange J) interaction almost by 
one order of magnitude. Qualitatively, a longer-range el- 
ph interaction results in a lighter mass of dressed carries 
(7 < 1) because the extended lattice deformation changes 

f dually in the process of tunnelling through the lattice 
[2l| . Also in the intermediate and strong-coupling 
regime the finite bandwidth is important, so that the ef- 
fective mass approximation cannot be applied. 

In this Letter we calculate isotope effects on the whole 
coherent band-structure and DOS of a two-dimensional 
lattice polaron with short- and long-range e-ph inter- 
actions by applying a continuous-time QMC algorithm. 
Unlike the strong-coupling limit J?] the isotope effect de- 
pends on the wave vector in the intermediate region of 
parameters. We also discuss the isotope effect on the 
electron spectral function including its incoherent part. 

Our model Hamiltonian describes any-radius e-ph in- 
teraction of an electron (or hole) on a square-lattice plane 
with linearly polarized vibrations of all ions of the 3D 
crystal plj . 
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FIG. 1: Top panels: small Holstein polaron band dispersions 
along the main directions of the two-dimensional Brillouin 
zone. Left: A = 1.1, right: A = 1.2. Filled symbols is lj = 
0.70 open symbols is = 0.66 Lower panels: the band 
structure isotope exponent for A = 1.1 (left) and A = 1.2 
(right). 



H = - ^[i(m - n)c]^Cn+./(m - n)^nc]^Cm]- 



q 



(1) 

Here Cm, and 5q are annihilation electron and phonon 
operators, respectively, i(m — n) is the hopping inte- 
gral, which is nonzero only for the nearest neighbors, 
t(a) = t (a is the primitive lattice vector), = 
(2ArA/tj)^i/2^^g»q n^^ _^ jj^ ^ jg ^j^g displacement op- 
erator at cite n, iV is the number of cells, and /(m — n) 
is a "force" applied to the electron at site m due to the 
ion displacement at cite n. 

In the strong-coupling, A ^ 1, and non-adiabatic 
limit Lo >, t one can apply the Lang-Firsov transfor- 
mation [23 to obtain the ground state energy Ep = 
(2Muj)~'^ 'Ylim ./^(™) and the coherent polaron band dis- 
persion Ck = i?k exp(— g'^) with the polaron narrowing 
exponent = [IMuj^Y^ Em [/'(m) - /(m)/(m -I- a)] 
and the "bare" band dispersion E\^. The band-structure 
isotope exponent a^, does not depend on the wave-vector 
k in this limit, because e-ph interactions do not change 
the band topology. 
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It is the same as a^. If the interaction is short-range, 
/(m — n) oc (5in,n (the Holstein model), then = Ep/u. 
Generally, one has — ^Ep/uj with the numerical coef- 
ficient 7 = (1 - Em /(m)/(m + a))/ p{n). 7 can 
be significantly less than 1 for a long-range (Frohlich) 
interaction po|. such as the unscreened interaction with 
c-axis vibrations of apex oxygen ions in the cuprates |2lj| : 



/(m — n) oc (|m — np + 1) 



-3/2 



(3) 



Here the distance between the in-plane and apex cites 
\/\m — np + 1 is measured in units of the in-plane lat- 
tice constant a, and the apex-plane distance is taken to 
be also a. Thus the strong-coupling isotope exponent, 
Eq. ((21) turns out to be numerically smaller for a longer- 
range e-ph interaction compared with a short-range in- 
teraction of the same strength Ep. 

QMC results (see below) show that this tendency also 
holds in the intermediate coupling regime, but the isotope 
exponent becomes a nontrivial function of the wave vec- 
tor: ah = Q!f,(k), because e-ph interactions substantially 
modify the band topology in this regime. We apply a 
continuous-time path-integral QMC algorithm developed 
by one of us 26j , which is free from any systematic finite- 
size and finite-time-step errors. The finite-temperature 
errors are exponentially small when the simulation tem- 
perature is smaller than the phonon frequency. The 
method allows for exact calculation of the whole polaron 
band dispersion on any-dimensional infinite lattice with 
any-range e-ph interaction using the many-body path in- 
tegral. 
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whre Co is the ground state energy. Configurations de- 
scribed by the electron position vector r and ion dis- 
placements ^, are obtained at imaginary time t = f3 — 
l/{kBT) from the configurations at r = by shifting 
along the lattice by vectors Ar. J VrV^^ integrates over 



3 




n 0.25 








0.2 








0.15 








. 0.1 








' 0.05 

i 















niiiTrTTTTTrTTrai| 







X M 

Momentum k 



2.5 p 

2 - 

1.5 ™ 

1 - 

0.5 - 

— 



X M r 

Momentum k 



Froehlicti A.=2.5 



Holstein X=1.2 




0.05 0.1 0.15 0.2 0.25 0.1 0.2 0.3 0.4 
{E-e„)/t (E-e„)/t 



FIG. 2: Top panels: small Frohlich polaron band dispersions 
along the main directions of the two-dimensional Brillouin 
zone. Left: A = 2.5, right: A = 3.0. Filled symbols is lj = 
0.70 1, open symbols is a; = 0.66 Lower panels: the band 
structure isotope exponent for A — 2.5 (left) and A = 3.0 
(right). 



single-electron paths with all possible shifts satisfying the 
twisted boundary conditions with the weight ^(r, ^). The 
statistics for any number of k points in the Brillouin zone 
are collected during a single QMC run [2^ . 

The isotope exponent ab(k) and DOS are presented in 
Figs. (1,2), and Fig. 3, respectively, for the small Holstein 
polaron (SHP) with the short-range interaction, and for 
the small Frohlich polaron (SFP) with the force given 
by Eq. The polaron spectra are calculated for two 
phonon frequencies, lo — 0.70 1 and uj — 0.66 whose 
difference corresponds to a substitution of O^^ by O^^ in 
cuprates. There is a significant change in the dispersion 
law (topology) of SHP, Fig.l, which is less significant for 
SFP, Fig. 2, rather than a simple band-narrowing (see also 
[26]). As a result, the isotope exponent 
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is k-dependent, Fig.l and Fig. 2 (lower panels). The 
strongest dispersion of ab is observed for SHP. Impor- 
tantly, the isotope effect is suppressed near the band edge 
in contrast with the k-independent strong-coupling iso- 
tope effect, Eq. It is less dispersive for SFP, espe- 
cially at larger A, where the polaron energy band is well 
described by the strong-coupling Lang-Firsov expression. 

Polaron DOS, N{E) = J^k 6{E-e]^), Fig.3, is less sen- 
sitive to topology, approximately scaling with the renor- 
malised bandwidth. It reveals a giant isotope effect near 
the van Hove singularity because the bandwidth changes 
significantly with to even in the intermediate-coupling 
regime. 

The coherent motion of small polarons leads to metal- 
lic conduction at low temperatures. Our results. Figs. 1-3, 



FIG. 3: Isotope effect on the polaron density of states, 
line is ci; = 0.70 i., dashed line is a; = 0.66 1. 
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show that there should be anomalous isotope effects on 
low-frequency kinetics and thermodynamics of polaronic 
conductors which depend on the position of the Fermi 
level in the polaron band. In fact, such effects have been 
observed in ferromagnetic oxides at low temperatures , 
and in cuprates 0, Q . To address ARPES isotope expo- 
nents one has to calculate the electron spectral func- 
tion A(k, E) taking into account phonon side-bands (i.e. 
off diagonal transitions) along with the coherent polaron 
motion (diagonal transitions). Using the 1/A expansion 
one obtains 12] 



A{k,E)=Y,[A\-\k,E) + A\+\k,E) 



(7) 



where 
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Here Z = exp(— E'p/w), ((k) = — fi, is the chemical 
potential, n{E) = [exp{l3E) + 1]'^, k;='= = k ± X;t=i Ir- 
and 7(q) is the Fourier transform of the force /(m) — 
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Clearly, Eq. Q is in the form of a perturbative multi- 
phonon expansion. Each contribution A^^\k, E) to the 
spectral function describes the transition from the ini- 
tial state k of the polaron band to the final state k^ 
with the emission (or absorption) of I phonons. Differ- 
ent from the conventional Migdal self-energy ^ the elec- 
tron spectral function comprises two different parts in the 
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strong-coupling limit. The first (/ = 0) k-dependent term 
arises from the coherent polaron tunnelling, ^co/i(k, E) — 



A[;~\'k,E) + A\j^'{k,E) = ZS{E - Ck) with a sup- 
pressed spectral weight Z <C f . The second incoherent 
part AincohO^i E) comprises all terms with ^ > 1. It de- 
scribes excitations accompanied by emission and absorp- 
tion of phonons. We note that its spectral density spreads 
over a wide energy range of about twice the polaron level 
shift Ep, which might be larger than the unrenormalised 
bandwidth 2zt in the rigid lattice. On the contrary, the 
coherent part shows a dispersion only in the energy win- 
dow of the order of the polaron bandwidth. It is im- 
portant that the incoherent background Aincoh (k, E) is 
dispersive (i.e. k-dependent) for the long-range interac- 
tion. Only in the Holstein model with the short-range 
dispersionless e-ph interaction 7(q) = const the incoher- 
ent part has no dispersion. 

Using Eq. {T)) one readily predicts an isotope effect on 
the coherent part dispersion and its spectral weight 
Z, and also on the incoherent background because Z , 
7(q), and lo all depend on M. While our predic- 
tion is qualitatively robust it is difficult to quantify the 
ARPES isotope effect in the intermediate region of pa- 
rameters. The spectral function, Eq. iQ), is applied in the 
strong-coupling limit A 3> 1. While the main sum rule 
(ii?A(k, = 1 is satisfied, the higher-momentum 

integrals dEEPA{k., E) with p > 0, calculated with 
Eq. lO differ from their exact values by an amount 
proportional to 1/A. The difference is due to partial "un- 



dressing" of high-energy excitations in the side-bands, 
which is beyond the lowest order 1/A expansion. The 
role of electronic correlations should be also addressed in 
connection with ARPES. While the results shown in Figs. 
1-3 describe band-structure isotope effects in slightly- 
doped conventional and Mott-Hubbard insulators with 
a few carriers, their spectral properties could be signifi- 
cantly modified by the polaron-polaron interactions 2^ , 
including the bipolaron formation [s^l at finite doping. 
On the experimental side, separation of the coherent and 
incoherent parts in ARPES remains rather controversial 

m 

In conclusion, we have calculated the isotope effect 
on the band structure of doped insulators employing a 
continuous-time quantum Monte-Carlo algorithm, and 
found its dispersion in the intermediate coupling regime 
which essentially depends on the radius and strength 
of the electron-phonon interaction. Using the strong- 
coupling electron spectral function we predicted an iso- 
tope effect on the weight and dispersion of the coherent 
part, and on the incoherent background. The exact iso- 
tope exponents. Fig. 1 and Fig. 2, could be instrumental 
in understanding of current and future experiments with 
isotope substituted oxides, and in assessing different an- 
alytical and semi-analytical approximations. 
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